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A dual acceptance criterion based on the sample mean and an extreme order is used in many inspection 
procedures. Computation of the acceptance probability for such a dual criterion is investigated. An ap- 
proximation and a lower bound to the acceptance probability are derived and are applicable to any con- 
tinuous distribution. In addition, the connection between this dual criterion and hypothesis testing of scale 
and location parameters is studied. In the case of the exponential distribution the exact evaluation of the 
acceptance probability yields the power of the test. 
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1, Introduction 

Suppose that a random sample of size n from a lot is measured with respect to a particular variable 
and that the acceptance or rejection of the lot depends upon whether or not the measurements satisfy 
certain criteria. "Lot" can refer to a group of individual items or to a specified amount of material 
which can be sampled randomly. 

There is widespread interest in sampling procedures that specify acceptance criteria involving- the 
sample mean and a proportion of defectives in the sample [1], [4], [S]. [9], [11] and [1 4]. 1 Such a sam- 
pling procedure might specify that the lot is to be accepted only if the sample mean is greater than a 
value n , say, and if no more than a specified percentage of the sample is less than a lower limit L. The 
purpose of a dual acceptance criterion is to ensure, for example, that the lot is at least a stated amount. 
^ u , of the specified variable on the average and that the number of so called "defectives" or items that 
violate the lower limit is controlled. Obviously, depending on the application, the acceptance criteria 
can be specified in the opposite direction; Le., the lot is to be accepted only if the sample mean is less 
than /J and at least a certain percentage of the sample is greater than an upper limit U. 

Specifically, let X 1 .*"* 1 X n be a random sample of n measurements, and let X (] > <■••< X (h| be the 
corresponding order statistics. It is assumed that the random variables X],"*,X n are independent and 
identically distributed (i.i.d.) with a probability density function fix), and that the X- have finite mean 
l* and variance o^. Let X be the sample mean and N^ be the number of defectives or measurements 
having values smaller than the specified (lowerl limit L. 

The sampling procedure to be considered is such that the lot is accepted whenever 

\X > pa and N L < k) (1.1) 

where ^ and k are specified in the sampling plan. 

In terms of the order statistics, 1 1 .1 1 is equivalent to the criterion 

|X^ J u andX[ Jj[+/( >i] 11.2) 

and the probability of accepting the lot is defined to be 

P n =P[X>^N L *Zk]. 11.3) 

The sampling procedure discussed above is a mixed variables-attributes acceptance criterion based 
on one sample. There are various ways of designing a mixed sampling plan. The type studied by Schill- 
ing and Dodge [19] is a double sampling procedure involving variables inspection in the first sample. 
If the variables inspection does not lead to acceptance 1 a second sample is taken and an attribute in- 
spection is conducted on the combined samples. In their work. Schilling and Dodge assume a normal 
distribution with unknown mean and known variance. 

We concentrate on a single sample plan where both the variables inspection as specified by the sam- 
ple mean and attributes inspection as specified by k, the number of allowable defectives, are con- 
ducted on the same sample. This causes difficulties in the computation of the acceptance probabilities 
because of the lack of independence of the sample mean and the order statistics. 

Investigations, of which we are aware, into the statistical properties of sampling procedures of this 
type assume a normal distribution with unknown mean and known variance, For instance In a com- 
pliance sampling application, Weed [21] simulates a two-stage procedure used in specifications for the 
thickness of paving material in which both stages involve a variable and an attribute inspection. Elder 
and Muse [8] develop a large sample approximation for the acceptance probability used in U.S. 
Department of Agriculture inspection procedures (1.3} and compare the approximation to an exact 
numerical procedure. 



Figures in brackets indicate literature reierences at the endoi this paper. 
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It is noted that the dual sampling criterion leads to an acceptance region for testing hypotheses con- 
cerning- the mean n and the probability of item defectiveness simultaneously. The probability of a 
defective is defined to hep = P[X *£ L). The acceptance region in (1.1) or (1.2) may be used for testing 
the null hypothesis 

Hq-.h = h andp=p 

versus the altern ati ves (1.4) 

i/ 1 :/i<M or p>p 

Through reparametrization, these hypotheses may be formulated in terms of the location and scale 
parameters. Evidently, this depends on the properties of the distribution under consideration. 
In the case of the normal distribution N\p, o 2 ), the probability of a defective is 



yher 



Thus, 



1 r z 

* b) = [ exp{-u 2 /2}du. 



a = [L - y.)/Q Hp\. (1.6) 

Consequently, u = p* andp =p* if and only if 

H = n* and a = a* = fL-jiVO -1 !/)*!. 
Accordingly, the hypothesis testing problem in (1.4) becomes that of testing 

i/ : ji = n and a = a 
versus (1.7) 

H V . M < ii* or a <^T^ } - 

Perusal of the literature turned up very few papers that are directly related to a joint test of the loca- 
tion and scale parameters. Eisenberger [7] develops an asymptotic joint test for the mean and variance 
of a normal distribution based on a quantile. Perng [18] develops a joint test for the location and scale 
parameters of an exponential distribution based on Fisher's method of combining two test statistics. 
Anderson [2] discusses the likelihood ratio test for simultaneously testing the mean and variance in 
multivariate normal distributions; both one-sample and /c-sample problems are considered. In a recent 
paper, Perlman [17] shows that the likelihood ratio testis unbiased. None of these papers discusses the 
computation of acceptance probabilities under alternative hypotheses. Also, unlike (1.7), the alter- 
natives in the quoted papers are rectangular regions. 

2. Scope of the Study 

It is our intention to investigate the acceptance probability of a dual sampling procedure from 
several aspects. The investigations are carried out for the normal distribution because of its im- 
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portance in acceptance sampling and for the exponential and Weibull distributions because of their 
application in modeling the life span distribution. 

First, in section 3, we derive a large sample approximation P a for the acceptance probability P„. 
This is achieved by deriving the asymptotic joint distribution of \/n{X —ft)/c and [lW L —np)/(np 
\\—p))^ as the sample size approaches infinity. This approximation method applies to any distribu- 
tion. We illustrate its use in the normal, Weibull. and exponential distributions. The results as given in 
sections 3.1, 3.2, and 3.3 are compared with a simulation study. 

In section 4 a lower bound P is established for P n that amounts to assuming the independence of the 
sample mean and the /c th order statistic, This lower bound for finite samples provides some informa- 
tion ou the accuracy of the approximation. We attempt to determine under what conditions the ap- 
proximation P a is a significant improvement over the lower bound. In this connection one notes that a 
large sample approximation P s is derived by normalizing the sample mean as \fn{X—y)/o and the 
number of defectives in the sample as W L — np)/(np\\ — pll- 1 '' 2 . If, instead, we convert N^ to an order 
statistic X|£| and consider X^ (or X^_ ki \ as an extreme statistic, the normalized sample 
mean \fn{X —\i) and X fle) (or equivalently X^^} are asymptotically independent (The proof is given in 
appendix B). This suggests that P_ serves as a possible approximation to P n when a is large and k is 
small. 

In other words, when comparing P^ and P one should keep in mind the relationship between k and 
n; namely, the ratio h/n. In the case of P s we have A^ /n ~* p and in the case of an extreme statistic we 
have kin -* as n -*■ °°. Clearly, one would expect that the lower bound P may be a reasonable ap- 
proximation when k is relatively small compared with n. This is indeed confirmed in our numerical 
study in section 4. The numerical studies show that P s is comparable to P for small kJn and superior 
to Pfor larger values of k/n. 

Finally, in section 5 the acceptance probabilities are approximated for the normal and Weibull 
distributions using a procedure proposed by Pearson and Hartley [16]. The exact acceptance prob- 
abilities curves are computed for the exponential distribution. 

3. Large Sample Approximation of the Joint Distribution of X andW^. 

3,1 Derivation 

LetX|,'",X n be a random sample from the lot with pdf ffx). Assume that X- has a finite mean \x and 
variance a'. 



Introducing indicator random variables /■, 
where 




13.1.1) 



and letting the probability that an item violates the lower specification limit L be 

p = P[X } - < L], (3.1.2) 

we can write the number of (unitl lower limit violations iV^in the sample as 

N L = Z /,. (3.1.31 

Note that N^ has a binominal distribution B{n,p\, and the event [A^ < k] is equivalent to the event 
[X^ +1 ) > L], In order to develop an approximation formula for the acceptance probability 

P n =P\X> M Q,N L <k\, 
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we consider random variables W n and Y n defined as 

W n = n i/2 (X-ii)/a 
and (3.1.4) 

Y n = {N L -np)/inp\l-p)) U2 . 

Let ( W n , Y n ) be a row vector. We prove the following result. 

THEOREM 3.1. As n— > «>, the random vector (W n , Y n ) converges in distribution to a bivariute nor- 
mal distribution with mean (0,0)' and covariance matrix 



(3.1.5) 



where 

c = E (X rM )I./o( P (l-p)) 1/2 . <3.1.6> 

PROOF: Let t j and t 2 be arbitrarily chosen but fixed real numbers. Form the linear combination of W n 
andK„, tl fF n +t 2 Y n . 

Direct computation and application of the central limit theorem give 

D 

t, W tl +t 2 Y ll — N (0, t] 2 + t 2 2 + r^afi) as n— «> 

It then follows from application of the Cramer-Wold device that 

- N\ ( ], I as n-«> 

• I \W ) 

where X is given in (3.1.5). 

Making use of the asymptotic distribution in Theorem 3.1, we note from (3.1.4) that 

X = n- l ' 2 oW n + u 
and 
N L = {nptt-p)) l/2 Y n +np. 

Thus the random vector {X , N[)' has asymptotically a bivariate normal distribution with mean and 
covariance matrix l~ given by 

/ o£ E(X r v)I f \ 

( ] andT = J (3.1.7) 

\n P J \£/X / -n»/ j npil-p) } 

respectively. 

For convenience in computation, write the acceptance probability P n as 

P n = P[X > Mo ] - P[X > Mo , N L > k] 

= P[X^ Mo ]-P[^ n >v^(M -M)/o,^ n >Ml-pH- 1/2 (A-«p)]. (3.1.8) 
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Making use of 13.1.7) and the continuity correction factor 0.5 for the random variable N L , we see that 
for sufficiently large n, P n may be approximated by 

CO OO CO 

P a = Jexp(-z 2 /2Mz- J" f g{x,y, Q )dxdy 13.1.9) 

\f2n & ah 

where 

a = VKlnff-nW a, {3.1. 10) 

b = (np(l-pir l/2 (/e+ 0.5 -r.pl, (3.1.11) 

g{x,y, Q ) = {2n)- i (l- e 2 r 1/2 e X p{-^ 2 4- r 2_2 eary )/2(l- e 2 )}, (3.1.12) 

and q is defined in 13.1.6). 

In order to compute the P [X > u , N^ < /c] using the approximation P s , we need to know the 
mean p and the variance o 2 of the distribution in question, the proportion defective p as defined in 
(3.1.2) and the correlation coefficient q as defined in (3.1.6). The computation of the bivariate normal 
term is described in more detail in Appendix A. 

3.2 Normal Distribution 

Assume that the sample comes from a normal distribution iV(ji,a"). 
The item defective probability from (3.1.2) ia 

p = P\X<L] = *{(L-^/o}, 13.2.1) 

where <t>{(L-(i)/cr} is the cdf of the ./V(0,1) given in(1.5). 

In order to compute the approximation P a given in (3.1.9), we need to compute the correlation coef- 
ficient given in (3.1.6K 

The expectation E[(X— pWfx^L]} ' s evaluated as 

E{(X-m)/ [X «L]1 = -~exp {-IL-h)V2o2J. 

Consequently the correlation coefficient is 

e = H2irp(l-p)r 1/2 exp Hl^)2/2o 2 }. 

In order to compare the approximation P a in (3.1.9) with an approximation developed by Elder and 
Muse [8], the lower limit L is chosen under the assumption that n = 0, a — 1, and according to the 
criterion 

P[N L ^k] = l-a, (3.2.2) 

where < a < 1 . 

Because A^ is B(n,p), the lower limiti is determined from 



; ?o (")"' 



p/(l-p>»-/'- 1-a, (3.2.3) 



where p = <t>(jL). 

Values of L as tabulated by Elder and Muse for a = 0.10, 0.05, and 0.01 are shown in table I. Once 
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TABLE I. Lower Limits used in Computation of Acceptance Probabilities for Normal Distribution 





k 


Lower Limit L 


n 


a = 0.10 


a=0.05 


o=0.01 


5 





2.036 


2.319 


2.877 




1 


1.215 


1.429 


1.843 




2 


0.685 


0.881 


1.250 


10 





2.309 


2.568 


3.089 




1 


1.602 


1.789 


2.157 




2 


1.196 


1.358 


1.670 


20 





2.559 


2.799 


3.289 




1 


1.928 


2.095 


2.428 




2 


1.586 


1.726 


2.001 


30 





2.696 


2.928 


3.402 




1 


2.100 


2.258 


2.574 




2 


1.783 


1.914 


2.172 



L is determined the correlation coefficient of X and N^ can be evaluated as 

e = -[2?rp(l-p)r 1/2 exp{-Z, 2 /2}. 



(3.2.4) 



The* Elder-Muse approximation along with their exact results are compared with the corresponding 
values of P a in table II where L is chosen such that a = 0.10. 

The comparison with the exact values derived in [8] shows that even for small sample size P a pro- 
vides an excellent approximation to the acceptance probability P n , and its effectiveness increases as k 
gets larger. When k = 0, the percent error in P a as compared to the exact results is approximately 3 
percent. For k = 1, it is about 1 percent and for k = 2, it is less than 1 percent. The percentage errors 
in both P a and the Elder-Muse approximation when p = are shown below. 





Percent Error in Approximations 




k 


= 


k = i 


k 


= 2 






Elder 




Elder 




Elder 


n 


P 


Muse 


P 

B 


Muse 


P 


Muse 


5 


3.3 


1.0 


1.0 


1.8 


0.6 


1.2 


10 


3.1 


0.6 


1.0 


1.0 


0.6 


1.2 


20 


3.0 


0.2 


0.8 


0.6 


0.6 


0.8 


30 


2.6 


0.2 


0.8 


0.8 


0.4 


0.6 



3.3 Weibull Distribution 

Assume that the sample X i ,...,X lt comes from a two parameter Weibull distribution W{\,8) with 
scale parameter k, shape parameter 8 and pdf 



f(x) = iO/k) (x/Al^'expf-U/A) 9 } for*>0,A>0, 6>0 
The mean and variance are 

p = Ar<i+i/e) 

and 

o 2 = A 2 {r<i+2/e>-[r(i+i/e)] 2 } 



(3.3.1) 

(3.3.2) 
(3.3.3) 
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respectively where !"(•) is the gamma function. 

For < 6 < 1, X has a decreasing failure rate (DFR) distribution; for 6 3= 1, X has an increasing 
failure rate (IFR) distribution. For further information see Johnson and Kotz [13]. 

In the case of the Weibull distribution, the proportion defective p is defined from (3.1.2) and (3.2.2) 
as 



p = lX<L]=l-exp{-(L/X) }. 



(3.3.4) 



The expectation 



EXI 



[X<L] 



B L 

~ / xU/A^'exp {- (x/X) 9 }dx 

A 



= Al{<L/A}<U/6} 



and I(c,d) is related to the incomplete l~ -function [1 2]. 

Combining (3.1.6), (3.3.4) and (3.3.5), we find that the correlation coefficient is 

Table II. Comparison of Approximation P fl with 
Elder -Mine Values for 
P\X»i** t N L <k] where P(N, <k)=0.90 
for Normal Distribution MO, 1 ) 



(3.3.5) 





k = 


*=1 


k = 2 










Elder 






Elder 






Elder 


n 


f 


Exact 


P a 


Muse 


Exact 


P a 


Muse 


Exact 


p » 


Muse 


5 


-.8 


0.035 


0.034 


0.032 


0.036 


0.036 


0.034 


0.037 


0.036 


0.037 




-.6 


0.087 


0.082 


0,085 


0.089 


0.088 


0.089 


0.089 


0.089 


0.091 




-.4 


0.180 


0.168 


0.181 


0.184 


0.181 


0.188 


0.185 


0.184 


0.189 




-.2 


0.318 


0.300 


0.323 


0.324 


0.320 


0.332 


0.326 


0.325 


0.333 




.0 


0.488 


0.472 


0.493 


0.496 


0.491 


0.505 


0.499 


0.496 


0.505 




.2 


0.659 


0.667 


0.663 


0.669 


0.672 


0.674 


0.671 


0.672 


0.674 




.4 


0.801 


0.814 


0.802 


0.811 


0.814 


0.811 


0.813 


0.814 


0.812 




.6 


0.899 


0.910 


0.899 


0.908 


0.910 


0.906 


0.909 


0.910 


0.906 




.8 


0.956 


0.963 


0.955 


0.962 


0.963 


0.959 


0.963 


0.963 


0.959 


10 


-.6 


0.027 


0.026 


0.026 


0.028 


0.028 


0.026 


0.028 


0.028 


0.027 




-.4 


0.098 


0.091 


0.097 


0.100 


0.098 


0.099 


0.101 


0.101 


0.101 




-.2 


0.252 


0.236 


0.253 


0.257 


0.252 


0.261 


0.260 


0-257 


0.264 




.0 


0.480 


0.465 


0.483 


0.490 


0.485 


0.495 


0.494 


0.491 


0.500 




.2 


0.713 


0.732 


0.714 


0.725 


0.735 


0.728 


0.731 


0.736 


0.733 




.4 


0.876 


0.897 


0.876 


0.888 


0.897 


0.887 


0.893 


0.897 


0.891 




.6 


0.956 


0.971 


0.956 


0.966 


0.971 


0.965 


0.969 


0.971 


0.967 




.8 


0.956 


0.994 


0.985 


0.992 


0.994 


0.991 


0.993 


0.994 


0.993 


20 


-.4 


0.034 


.0.032 


0.034 


0.035 


0.034 


0.034 


0.036 


0.035 


0.034 




-.2 


0.174 


0.162 


0.174 


0.178 


0.173 


0.178 


0.180 


0.177 


0.181 




.0 


0.474 


0.460 


0.475 


0.483 


0.479 


0.486 


0.488 


0.485 


0.492 




.2 


0.781 


0.811 


0.781 


0.795 


0.814 


0.795 


0.802 


0.814 


0.802 




.4 


0.937 


0.963 


0.937 


0.950 


0.963 


0.950 


0.956 


0.963 


0.955 




.6 


0.981 


0.996 


0.981 


0.991 


0.996 


0.991 


0.994 


0.996 


0.993 


30 


-.4 


0.013 


0.012 


0.013 


0.013 


0.013 


0.013 


0.014 


0.014 


0.013 




-.2 


0.127 


0.118 


0.127 


0.130 


0.126 


0.129 


0.131 


0.129 


0.131 




.0 


0.470 


0.458 


0.471 


0.479 


0.476 


0.480 


0.484 


0.482 


0.487 




.2 


0.824 


0.861 


0.824 


0.839 


0.863 


0.839 


0.847 


0.863 


0.846 




.4 


0.958 


0.986 


0.958 


0.972 


0.986 


0.972 


0.978 


0.986 


0.977 




.6 


0.985 


0.999 


0.985 


0.995 


0.999 


0.995 


0.997 


0.999 


0.997 



% = o 
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e = [U{iL/» 9 , i/e}-w}/o{ P a-p)) i/2 

where \i and o are defined by (3.3.2), and (3.3.3) respectively. 
The limits of integration for the approximation (3.1.9) are 



(3.3.6) 



a = 



n 1/2 [^ -Ar( 1+1/0)] 

A{ni+2/e)-[ni+i/e)] 2 } 1/2 



(3.3.7) 



and b as defined in (3.1.11). 

As is the case in the normal distribution, the lower limit L is determined according to (3.2.3) and 
(3.3.4) for specified values of k and a. 

Explicitly 



L=A[-log e (l-p)] 1/9 . 



(3.3.8) 



The proportion defective p is tabulated in table III for a = 0.10, 0.05 and 0.01, n = 5,10, 20, 30 and 
A =0, 1,2,3, Corresponding lower limits!/ where A = 1 are shown in table IV. 



TABLE III. Proportion Defectives p used in Computation of 
Acceptance Probabilities 







Proportion Defeetivep 


n 


k 


o=0.10 


o=0.05 


o=0.01 


5 





0.0208 


0.0102 


0.00200 




1 


.112 


.0765 


.0330 




2 


.247 


.1890 


.106 




3 


.416 


.3425 


.222 


10 





0.0105 


0.00511 


0.00100 




1 


.0543 


.0365 


.0155 




2 


.1153 


.0870 


.0473 




3 


.1875 


,1500 


.0930 


20 





0.00525 


0.00256 


0.000500 




1 


.0269 


.0180 


.00759 




2 


.0564 


.0422 


.0227 




3 


.0902 


.0713 


.0435 


30 





0.00350 


0.00171 


0,000335 




1 


.0178 


.0120 


.00500 




2 


.0373 


.0278 


.0149 




3 


.0594 


.0468 


.0285 



The approximation P fl is compared to a simulation study where the acceptance probability was com- 
puted from 5,000 random samples. Simulation for the Weibull distribution was done by generating in- 
dependent uniform random deviates U i using a congruential random number generator and making 
the transformation 



X^At-log^-) 



1/9 



TheX,- are independent W{k,Q) r.v.s with pdf as shown in (3.3.1). 

Values of P a and simulated acceptance probabilities are tabulated in table V for Weibull distribu- 
tion mi,8\ for B = 1,2, 3.5. 
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The accuracy of the approximation P a as gauged by the simulation results is dependent on several 
factors; i.e., namely, the value of the shape parameter 9; a, the probability that the sample will contain 
more than the allowable number of defectives; n, the size of the sample; and ft, the number of 
allowable defectives or number of measurements less than the Lower limit L . 

The worst accuracy is for a Weibull distribution with 9 = 1 where a is small, a = 0.01, and n is 
small, n = 5. The error is 9 percent for this case but drops to 2 percent when the sample size is in- 

T.\BLE IV. Lower-Limits Used in Compulation of Acceptance Probabilities for Weibull Distribution 





n 


k 


Lower Limit/. 




» = 0.10 


»=0-05 


.7=0,01 


0=1 


5 





0.0210 


0.0103 


0.0020 




5 


1 


.1188 


.079b 


.0336 




5 


2 


.2837 


.2095 


.1120 




5 


3 


.5379 


.4193 


.2510 




10 





.0106 


.005 1 


.0010 




10 


1 


.0560 


.0372 


.0156 




10 


2 


.1227 


.0910 


.0487 




10 


3 


.2076 


.1625 


.0976 




20 





.0053 


.0026 


.0005 




20 


1 


.0273 


.0182 


.0076 




20 


2 


.0581 


.043 L 


.0230 




20 


3 


.0945 


.0:40 


.0445 




30 





.0035 


.0017 


.0003 




30 


1 


.0180 


.0121 


.0050 




30 


2 


.0380 


.0282 


.0150 




30 


3 


.0612 


.0479 


.0289 


9=2 


5 





.1450 


.1013 


.0447 




5 


1 


.3446 


.2821 


.1832 




5 


2 


,5326 


.4577 


.3347 




5 


3 


.7334 


.6475 


.5010 




10 





.1027 


.0716 


.0316 




10 


1 


.2367 


.1928 


.1250 




10 


2 


.3503 


.3017 


.2206 




10 


3 


.4557 


.4031 


.3124 




20 





.0726 


.0506 


.0224 




20 


1 


.1631 


.1348 


.0873 




20 


2 


,2409 


.2076 


.1515 




20 


3 


.3075 


.2720 


.2109 




30 





.0592 


.0414 


.0183 




30 


1 


.1340 


.1099 


.0708 




30 


2 


.1950 


.1679 


.1225 




30 


3 


.2475 


.2189 


.1700 


6 = 3.5 


5 





.3317 


.2702 


.1694 




5 


1 


.5441 


.4852 


.3791 




5 


2 


.6977 


.6398 


.5351 




5 


3 


.8376 


.7301 


.6737 




10 





.2724 


.2216 


.1390 




10 


1 


.4390 


.3904 


.3047 




10 


2 


.5492 


.5042 


.4216 




10 


3 


.6382 


.595 


.5144 




20 





.2233 


.1818 


.1140 




20 


1 


.3573 


.3181 


.2482 




20 


2 


.4434 


.4073 


.3402 




20 


3 


.5097 


.4752 


.4109 




30 





,1968 


.1620 


.1017 




30 


1 


.3171 


.2831 


.2202 




30 


2 


.3929 


.3607 


.3013 




30 


3 


.4502 


.4198 


.3633 
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Table VT Comparison of Approximation P with Simulation {or 

P(X**H a +. N L < k) where P(N L £k)=l-ir 

for Weibull Distribution Wil.8). 





n 


k 


Probability of Acceptance 




o = 


= 0.10 


o = 0.05 


o=0.01 




P 


Simul 


P Simul 

a 


P Simul 

a 


e=i 


5 





0.645 


0.634 


0.698 0.663 


0.712 0.672 




5 


1 


.668 


.652 


.698 .673 


.712 


672 




5 


2 


.678 


.664 


,699 .680 


.711 


674 




5 


3 


.688 


.674 


.702 .684 


.711 


675 




10 





.705 


.707 


.768 .755 


.785 


773 




10 


1 


.726 


.724 


.768 .763 


.785 


776 




10 


2 


.733 


.734 


.766 .768 


.785 


778 




10 


3 


.738 


.743 


.767 .772 


.784 


779 




20 





.775 


.795 


.848 .828 


.868 


875 




20 


] 


.795 


.801 


.846 .836 


.868 


877 




20 


2 


.798 


.807 


.844 .840 


.867 


879 




20 


3 


.801 


.810 


.843 .841 


.867 


878 




30 





.815 


.826 


.893 .872 


.914 


914 




30 


1 


.835 


.837 


.890 .879 


.914 


918 




30 


2 


.83 7 


.838 


.887 .878 


,914 


915 




30 


3 


.838 


.841 


.886 .881 


.913 


917 


6=2 


5 





.681 


.703 


.720 .731 


.744 


738 




5 


1 


.709 


.723 


.729 .734 


.744 


739 




5 


2 


.721 


.736 


.733 .737 


.744 


740 




5 


3 


.729 


,740 


.736 ,739 


.744 


740 




10 





.743 


.764 


.788 ,807 


,824 


822 




10 


1 


.768 


.780 


.803 .808 


.824 


825 




10 


2 


,777 


,794 


.808 ,808 


.823 


826 




10 


3 


.784 


.800 


.812 ,810 


,823 


827 




20 





.810 


.827 


.865 .885 


.906 


903 




20 


1 


.832 


.839 


.876 ,884 


.906 


905 




20 


2 


.837 


.848 


.884 .882 


.905 


907 




20 


3 


.840 


.851 


.884 .882 


.905 


907 




30 





.844 


,860 


.897 ,924 


.946 


947 




30 


1 


,865 


,866 


.903 .922 


.946 


950 




30 


2 


.367 


.867 


.909 .919 


.946 


949 




30 


3 


.869 


.872 


.914 .918 


.945 


952 


e=3.5 


5 





.801 


.829 


.864 .859 


.880 


872 




5 


1 


.830 


.844 


.865 .873 


.880 


875 




5 


2 


.839 


.855 


.868 .875 


.880 


876 




5 


3 


.843 


.853 


.869 .878 


.879 


877 




10 





.853 


.864 


.931 .911 


.951 


942 




10 


I 


.877 


.868 


.930 .915 


.951 


947 




10 


2 


.882 


.884 


.928 .922 


.951 


946 




10 


3 


.885 


.889 


.928 .925 


.951 


950 




20 





.882 


.890 


.968 .938 


.990 


983 




20 


1 


.902 


.894 


.964 .946 


.990 


984 




20 


2 


.903 


.899 


.960 .946 


.990 


987 




20 


3 


.903 


.894 


,958 .947 


.989 


986 




30 





.887 


.893 


,974 .938 


.998 


984 




30 


1 


.908 


.896 


.970 .947 


.998 


989 




30 


2 


,907 


.894 


.966 .948 


.997 


990 




30 


3 


.907 


.898 


.964 .956 


.996 


990 



*H ti = 0.75 



495 



creased to n = 1 0. For other Weibull distributions and combinations of a and n, the worst accuracies 
occur when k = 0, and in this case the errors are as large as 6 percent for n = 5 and 4 percent for n = 
30. However, the approximation P a works very well when h > 0. The disagreement between P a and 
the simulation is less than 1 percent for a large proportion of the points when k > 0. 

3.4 Exponential Distribution 

Assume that the sample X ]t ...,X H comes from an exponential distribution EiX,p) with location 
parameter p and scale parameter A and pdf 

f(x) = U/MexpHx-plA) x>M>0 &A.1) 

The mean and variance of X are given by u = A + /3 and o 2 = A~ respectively. 
We have 

p = l-ex P HL-/3)/A) (3.4.2) 

and 

EXI [X^L] = ^P~ tt-p)\L-p) + pp. (3.4.3) 

Combining (3.4.2) and (3.4.3), we get 

o=-(l-p) 1/2 (L-/3)/Ap 1/2 . (3.4.4) 

Using values for the proportion defective p that are given in table III, the corresponding limits L as 
determined by 

L=/J-Alog(l-p) (3.4.5) 

are found in table VI for /3 = and A = 0.5, 1 , 2. 

The values a and b appearing in the approximation P a (3.1.9) are given by 

a = n l/2 r l (m o -A-0> (3.4.6) 

b = inpil-p))- i/2 (h+ 0.5 ~np) 

and o is defined by (3.4.4.) 

Values of P u and simulated acceptance probabilities are tabulated in table VII for the exponential 
distribution E(X,0) for A = 0.5,1,2. 

The accuracy of the approximation P a is more dependent on n, the sample size and less dependent 
on A, the number of allowable defectives for the exponential distribution than for Weibull distribu- 
tions. The worst accuracy-is for an exponential distribution with A = 1, where k = and n = 5. The 
disagreement with the simulation in this case is 7 percent, dropping to 1 percent when the sample size 
is increased ton = 10. In general, the accuracies are not dependent upon the parameter A but are 
somewhat dependent upon the way in which the lower limit L is chosen, and the accuracies tend to 
worsen as the probability of the sample containing more than the allowable number of defectives in- 
creases. Accuracies of about 2 percent are characteristic of the results over all values of A. 

4. A Lower Bound for the Acceptance Probability 

A lower bound for the acceptance probability is provided by the following lemma. 
LEMMA 4.1: Let Xj,...,X n be Ltd random variables from a continuous distribution. Let X be the 
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sample mean and X( r ) be the r th smallest order statistic of X lv ..,X n 
real numbers a, b and positive integer r , 1 < r <n, 

PfX > a, X (r) 5= b] 3> P[X 3> a] P[X (r) > b] 

P[X < a, X (r) < b] > P[X < a] P[X W < b]. 



Then for arbitrarily fixed 

(4.1) 
(4.2) 



The lemma is an easy eonsequenee of a general theorem (Esary, Proschan, and Walkup [10]). For easy 
reference, we quote the theorem below, as well as the definition of "associatedness." Random 

TABLE VI. Lower Limits used in Computation of Acceptance Probabilities for Exponential Distribution 











Lower Limit L 




n 


k 


o=0.10 


o = 0.05 


o=0.01 


A=0.S 


5 





0.0105 


0.0051 


o.ooio 




5 


1 


.0594 


.0398 


.0168 




5 


2 


.1418 


.1047 


.0560 




5 


3 


.2689 


.2097 


.1255 




10 





.0053 


.0026 


.0005 




10 


1 


.0280 


.0186 


.0078 




10 


2 


.0614 


.0455 


.0243 




10 


3 


.1038 


.0813 


.0488 




20 





.0026 


.0013 


.0003 




20 


1 


.0136 


,0091 


.0038 




20 


2 


.0290 


.0216 


.0115 




20 


3 


.0473 


.0370 


.0222 




30 





.0018 


.0009 


.0002 




30 


1 


.0090 


.0060 


.0025 




30 


2 


.0190 


.0141 


.0075 




30 


3 


.0306 


.0240 


.0145 


A=1.0 


5 





.0210 


.0103 


.0020 




5 


I 


.1188 


.0796 


.0336 




5 


2 


.2837 


.2095 


.1120 




5 


3 


.5379 


.4193 


.2510 




10 





.01 06 


.0051 


.0010 




10 


1 


.0560 


.0372 


.0156 




10 


2 


.1227 


.0910 


.0487 




JO 


3 


.2076 


.1625 


.0976 




20 


(1 


.0053 


.0026 


.0005 




20 


1 


.0273 


.0182 


.0076 




20 


2 


.0581 


.0431 


.0230 




20 


3 


.0945 


.0740 


,0445 




30 





.0035 


.0017 


.0003 




30 


1 


.0180 


.0121 


.0050 




30 


2 


.0380 


.0282 


.0150 




30 


3 


.0612 


.0479 


.0289 


k=2 


5 





.0420 


.0205 


.0040 




5 


1 


.2376 


.1592 


.0671 




5 


2 


.5674 


.4190 


.2241 




5 


3 


1.0757 


.8386 


.5021 




10 





.0211 


.0102 


.0020 




10 


1 


.1121 


.0744 


.0312 




10 


2 


.2455 


.1820 


.0973 




10 


3 


.4153 


.3250 


.1952 




20 





.0105 


.0051 


.0010 




20 


1 


.0545 


.0363 


.0152 




20 


2 


.1161 


.0862 


.0459 




20 


3 


.1891 


.1479 


.0889 




30 





.0070 


.0034 


.0007 




30 


1 


.0359 


.0241 


.0100 




30 


2 


.0760 


.0564 


.0300 




30 


3 


.1225 


.0959 


.0578 
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TABLE VII. Comparison of Approximation P a with Simulation for P[X >fi u , 

for Exponential Distribution ECKO) 



"l< 



kl where P[N L <k]= l-a 





Probability of Acceptance 




o=0J0 


cr = 0.05 


cr=0.01 


n 


k 


P 

a 


Sim ill 


P Simul 

a 


P 

a 


Simul 


A= 0.5, jj = 0.25 














5 





0.781 


0.825 


0.850 0.853 


0.868 


0,893 


5 


] 


.804 


.846 


.847 .871 


.868 


.896 


5 


2 


.811 


.857 


.846 .874 


.867 


.897 


5 


3 


.819 


.869 


.847 .885 


.866 


.899 


10 





.842 


.875 


.922 .929 


.943 


.957 


10 


1 


.862 


.889 


.919 .927 


.943 


.958 


10 


2 


.865 


.893 


.915 .933 


.942 


.958 


10 


3 


.807 


.897 


.913 .935 


.941 


.956 


20 





.878 


.902 


.964 .948 


.987 


.985 


20 


1 


.898 


.905 


.960 .954 


.987 


.984 


20 


2 


.898 


.907 


,956 .955 


.986 


.984 


20 


3 


.898 


.899 


.953 .948 


.985 


.984 


30 





.886 


.901 


.973 .952 


.997 


.987 


30 


1 


.906 


.904 


.969 .951 


.997 


.990 


30 


2 


.906 


.903 


.965 .948 


.996 


.990 


30 


3 


.906 


.903 


.962 .950 


.995 


.989 


A= 1.0, jj = 0.75 














5 





.645 


.632 


.698 .651 


.712 


.686 


5 


1 


.668 


.657 


.698 .660 


.712 


.689 


5 


2 


.673 


.665 


.699 .667 


.711 


.690 


5 


3 


.688 


.676 


.702 .674 


.711 


.691 


10 





.705 


.712 


.768 .754 


.785 


.776 


10 


I 


.726 


.726 


.768 .761 


.785 


.777 


10 


2 


.733 


.734 


.767 .766 


.785 


.777 


10 


3 


.738 


.745 


.767 .769 


.784 


.778 


20 





.775 


.789 


.848 .846 


.868 


,873 


20 


1 


.795 


.800 


.846 .849 


,868 


,876 


20 


2 


.798 


.804 


.844 .849 


,867 


,875 


20 


3 


.801 


.806 


,843 .851 


.867 


.873 


30 





.815 


.839 


.893 .890 


.915 


.920 


30 


1 


.835 


.846 


.890 .887 


.914 


.921 


30 


2 


.837 


.842 


.887 .884 


.914 


.922 


30 


3 


.838 


.846 


.886 ,891 


.913 


.921 


X= 2.0, jj = 0.75 














5 





.825 


.872 


.899 .915 


.919 


.950 


5 


1 


.846 


.883 


.895 .927 


.919 


.954 


5 


2 


.851 


.891 


.893 .929 


.918 


.953 


5 


3 


.856 


.894 


.892 .935 


.916 


.955 


10 





.870 


.894 


.954 .944 


.976 


.987 


10 


1 


.889 


.889 


.950 .944 


.976 


.987 


10 


2 


.891 


.898 


.945 .948 


.975 


.986 


10 


3 


.891 


.894 


.942 .945 


.974 


.985 


20 





.887 


.902 


.974 .952 


.997 


.991 


20 


1 


.906 


.903 


.970 .956 


.997 


.989 


20 


2 


.906 


.905 


.965 .956 


.996 


.991 


20 


3 


.906 


.904 


.962 .955 


.995 


.992 


30 





.889 


.907 


.976 .956 


1.000 


.990 


30 


1 


.909 


.906 


.972 .949 


.999 


.991 


30 


2 


.908 


.902 


.968 .948 


.999 


.991 


30 


3 


.908 


.901 


.965 .952 


.998 


.990 



498 



variables X ],..., X n are said to be associated if 

Cov(f(T),g(T))>0 

for all non-decreasing functions / and g in each X- for which Ef(T), EgfT), Ef(T)glT) exist and T 
denotes {X : ,...,X„}. 

THEOREM 4.1. Let T = {X 1 ,...,X n } be associated, Sj = f;(T) and fj be nondecreasing for i=l„..,k. 
Then 

k 

PfS^sj S k <s k ]> n P[Si<5,.] (4.3) 

'-1 

k 

P[S]>Sj S k >s k ]> n P[5,->5,.] 14.4) 



i=l 



for all s^,...,sy. 



PROOF OF LEMMA 4.1: In our case the X ; 'a are statistically independent and hence associated. Let 
5; = X and S 2 = X( r )- Clearly, 5j and S 2 are non-decreasing functions in each of the X ( .'s; hence (4. 1) 
and (4.2) hold. Moreover, Cov(Sj, S 2 ) = Cov (X , Xjy,) 3= 0. This completes the proof. 

From Lemma 4. 1 , we have a lower bound P to the acceptance probability 

P=P[X >a)P[X (k+1) >L}<P n =P[X >a,X ik+n >L], (4.5) 

where k + 1 corresponds to r. 

The r.v. X^ +l ^ can be transformed to a r.v. Z with Beta distribution with parameters n— k and A+l. 
Thus 

rin+1) \~p 

P[X (k+I) >L] = P[Z<\~F x iL)} = / *»-*-! (l-«)*<fe. (4.6) 

r(A+l)r|n-A) 

The lower bound P in (4.5) can be computed using the marginal distribution of the sample mean and 
the Beta distribution. 

Because the computation of the lower bound P is much easier than the computation of the accep- 
tance probability P n it would be an immense simplification if the lower bound could serve as an ap- 
proximation for P„ . 

Therefore, it is of practical importance to determine the samplesize n and values of k that are 
necessary in order that the lower bound be an acceptable approximation for P n , In other words, it is of 
interest to know the smallest value of n and the range of k values which makes the independence of X 
and X( )c+1 ) acceptable. 

5. Comparison of the Exact Probability of Acceptance with the 
Approximation and the Lower Bound 

5.1 Acceptance Probability Curves 

The acceptance probabilities computed using either simulation or numerical integration along with 
the corresponding lower bound P and the approximation P a are plotted as a function of one parameter 
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of the distribution in question. This provides a comparison of the relative accuracy of P a to P as a 
technique for approximating P n . The curves are varied over n and A: in order to examine the effect of 
sample size and number of allowable defectives A: on P n ,P a and P. 

5.2 Normal Distribution 

Assuming that X t ,..., X n are i.i.d. N(\i,\ ), the acceptance probability 

P u = P /l \X>l' ,X (k+1 ,>L] 

for L chosen according to (3.2.3) and 1*0=0 was computed using a technique for simulating random 
normal deviates due to Box and Muller [3]. The resulting acceptance probabilities as a function of n 
are shown as the solid line in figures 1 -4. 

The corresponding lower bound P was computed from (4.5) and the approximation P a was com- 
puted for (3.1.9). 

The relationships among the probability of acceptance P n , its approximation P a , and its lower 
bound P as a function of sample size n and allowable number of defectives k is depleted in figures 1-4 
for samples of size n = 10 and n = 30. The following convention is used for all figures; namely, P n is 
shown as a solid line; P a is shown as a heavy dashed line; and .Pis shown as a lighter dotted line. 

From figure 1 it is obvious that when A: = and n is small, P a is a better approximation to the ac- 
ceptance probability than the lower bound as long as \i < 0.25. As n increases the superiority of P to P_ 
increases as A: is allowed to become larger. For example, when k = 3 as in figure 4, the lower bound 
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Approximation f^ dr 


- 


>f'Lower bound P 




J? 







-e.s 
(a) 




FIGURE 1 . Acceptance probabilities when the number of allowable defectives k = and n observations are drawn 
from the normal distribution JV(fj, 1 ). 
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Figure 2. Acceptance probabilities when the number of allowable defectives k = 1 and n observations are drawn 
from the normal distribution JVijJ, 1). 
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Figure 3. Acceptance probabilities when the number of allowable defectives k = 2 and n observations are drawn 
from the normal distribution N[fi, 1 ). 





Figure 4. Acceptance probabilities when the number of allowable defectives k = 3 and n observations are drawn 
from the normal distribution Niyi, 1 1. 

does not give a satisfactory approximation for the smaller sample size, and P fl is clearly preferable. 
Even for n = 30, P a is at least as accurate as P over the entire range of u. 

5.3 Weibull Distribution 

Assuming that X 1 ,...,X rl are i.i.d. fF{l,6), and that ^ = 0.75 and that L is chosen according to 
(3.3.8) with 6=1, the acceptance probability was computed by simulation and is shown as the solid 
line in figures 5-8. The corresponding lower bound P was also computed using simulation and is 
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FIGURE 5. Acceptance probabilities when the number of allowable defectives k = and n observations are drawn 
from a Weibull distribution W{\, 6). 
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FIGURE 6. Acceptance probabilities when ihe Dumber of allowable defectives k = 1 and n observations are drawn 
from a Weibull distribution W{\, 8). 
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FIGURE 7. Acceptance probabilities when the number of allowable defectives k = 2 and n observations are drawn 
from a Weibull distribution JT(1, fl). 
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FIGURE 8. Acceptance probabilities when the number of allowable defectives k = 3 and n observations are drawn 
from a Weibull distribution W{1, 9). 

represented by the dotted line in the same figures. The approximation P a is shown by the heavy 
dashed line in the figures. 

The figures show that P a is not a particularly good approximation to P n when k = 0, and one would 
do much better using the lower bound P. However, P a shows the same characteristic for the Weibull 
distribution as for the normal distribution; namely, that as fc/n increases the accuracy of the approx- 
imation increases. For n = 10 and k = 3, P a is superior to P} for n = 30, Pis indistinguishable from 
the simulated acceptance probability. 
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5.4 Exponential Distribution 

5.4.1 Comparison with a UMP Test 

As discussed in section 1, we may view the problem of finding an optimal sampling procedure as a 
hypothesis testing problem formulated in (1.4). In general there exists no uniformly most powerful 
(UMP) test for (1.4). However, it is interesting to note that in the exponential distribution the dual ac- 
ceptance criterion for k = corresponds to a test which is UMP for a subset of alternatives specified 
in (1.4). Specifically, suppose the sample comes from the exponential pdf given in (3.4.1). 

The UMP acceptance region for testing 

H : A = A* and (i = 0* 
versus 

#! : 0< A < A* and 0< /3 < /?* 
is given by 

[X >K,,X m >/}*]. (5.4.1) 

This testing problem is equivalent to testing 

Hq: A = A* andp = p* 
versus 

J ff 1 :0<A<A* andp > 1 - (1 - p *)iVA 
where 

p* = l-exp{-(L-/»*)/A*} 
or 

|3* = L + A*Iog(l-p*). 
Under if(j, Px*,p* I^U) ^ 0*1 = ^' an ^ Mo * s determined by the equation 

P i , i p,|X>Mol=l-<». ' 5 - 4 -2l 

where a is a predetermined level of significance (Lehmann [15]). 

If we set L = fi* and k = 0, the test specified by (5.4.1) clearly is the same test specified by (1.3), 
and the acceptance probability 

can be computed either by the approximation shown in section 3.4 or by numerical integration using 
an exact formula for the distribution of X and N^ as shown in the next section, 

5.4.2 Exact Distribution of X andNj 

The joint distribution of X and N L can-be obtained from the order statistics. 
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Let 



We have the pdf of Z< i ), 



and for i > 2, Z ; has a pdf 



Z, — nX (1) 
Z,. = («-,-+ DIX,,.,-*^,). 

gjbi) = A" 1 exp {-{zi-nffl/li},zi>p (5.4.4* 

*,.(*,•> = A" 1 exp (^ ; /A), £i .>0. 



To compute the acceptance probability P for an arbitrary k, we make use of the fact that the Z's 
are independent r. v. 's, and that 

n n n 

X X-, = Z Z ; = X X ( ;k and proceed as follows: 

P n =P Xfi [X>^N L <k\ 

n k+\ , fA+1 | 

= / »•//» [JZ,.^ n Mo , Z Z,7(n-i+ 1 ) > L |*i ,-, **+ 1] 1 n ft U,.) } <fe, ~dz k+ j 

re ft+1 ) jfc+1 ) 

= ; 7 ; p [ *+ 2 Z; * " Mo " f * fl i i g '■ (z '■ , ' ^i"* dz *+i (5 * 4 - 5) 

where A — {fej,"', ^ + j): £ z r -/(n-(+ 1) > L and nfio ~~ 2z ( > 0|. 

The expression in (5.4.5) is the exact probability of acceptance, P n . 

When k = 0, the computation of P n reduces to 

P n = T ° P [ I Z { > n Mo - Z] ] g^zjldzj+f gfatej. (5.4.6) 

ni ^ n)J 

n 

Note that the sum Y = 2. Z i has a gamma density. 

/(y) = (1/A)"" 1 y~ 2 exp(-y/A) . (5.4.7) 

r(n-l) 

Substituting (5.4.4) and (5.4.7) in (5.4.6) we obtain 

P„ = I J jV" v n - 2 exp{~{z l -nP)/X}dvdz 1 + exp {-n| Mo -/J)/A} (5.4.8) 

Hn-1) a c 



where 



The lower bound for P„ is 



a = nL 

b = nfio 

c = (/ifi -zj)/A 



00 QO 

- ~ Jfx)dx JgMdx 



504 



(5.4.9) 



where fix) is the pdf of the x 2 (2n) and gipc) is the pdf of the x 2 '2). 

5.4.3 Acceptance probabilities 

If we assume that X j, "*,X n are i.i.d. £U,0), the acceptance probability for/c = 

P n = P K {X > „ , N L < 0] = P\X > mo- *iu > IA 

is computed from (5.4.8) using a numerical integration technique that takes advantage of the fact that 
the inner integral is an incomplete T-function. Note that /Jq is determined from x 2 (2l according to 
{5.4.2), and L is determined according to (3.4.5). The acceptance probability P n is shown as the solid 
line in figure 9. 





FIGURE 9. Aceeptanee probabilities when the number of allowable defectives k =- and n observations are drawn 
from an exponential distribution E{k, 0). 



The acceptance probabilities fork = 1,2,3, for/j^ = 0.75 and L chosen according to (3.4.5) were com- 
puted by simulation as were the corresponding values of P. The approximation P a was computed from 
(3.4.6). Results are shown in figures 10-12. 

The graphs show that P is a better approximation to P n than the lower bound Pfor small sample 
size where the superiority of P a over P_ increases as k increases. For large sample size, say n = 30, the 
two methods give almost identical approximations to P n . 



Values of loused in Computation of 

Acceptance Probabilities for TJMP Test for 

Exponential Distribution 

Values of n„ 



n 


o=0.10 


a=0.05 


a=0.01 


5 


0.48652 


0.39403 


0.25582 


10 


0.62213 


0.64254 


0.41302 


20 


0.77626 


0.66273 


0.55411 


30 


0.77431 


0.71998 


0.62475 
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FIGURE 1 1 , Acceptance probabilities when the number of allowable defectives k = 2 and n observations are drawn 
from an exponential distribution EU. 0). 
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FlGUBE 12. Acceptance probabilities when the number of allowable detectives/: = 3 and n observation 8 are drawn 
from an exponential distribution EiX, 01. 

6. Synopsis 

The problem of computing the acceptance probability P n has been addressed by an approximation 
P a that relies on the asymptotic Joint distribution of the sample mean and number of defectives in the 
sample. P B has the advantage that it ia applicable to any continuous distribution. It is computed using 
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a iV(0,l) cdf and a bivariate normal cdf which in turn can be reduced to a single variable integration. 

The approximation P a compares very favorably with another published approximation for the nor- 
mal distribution and with a lower bound P. Graphs of the acceptance probability as a function of one 
parameter of the distribution are used to compare the relative accuracies of P a and P. The graphs 
show that for the normal distribution P a and P have comparable accuracies with k = 0. As k/n in- 
creases, P a quickly becomes superior to P, and even for large n and k > P a is superior. In other 
words, the best results for the normal distribution are obtained with P when k = and with P a for all 
other values of k . 

In the case of Weibull distribution P is superior for k = 0. As k/n increases, P a gains in accuracy, 
and for large n, P continues to have an edge over P B . The difficulty in computing P for the Weibull 
distribution may make it desirable to use P B for all applications. 

In the case of exponential distribution, the exact joint distribution of the sample mean and number 
of defectives in the sample has been derived for k = 0. The computation of the acceptance probability 
P n in this case involves a two-variable integration. Graphs of the acceptance probabilities show that 
the lower limit P gives a consistently good approximation to the acceptance probability. The approx- 
imation P a and the lower limit P have also been computed for the exponential distribution for 1 < k < 
3. The graphs for these tests show that P is comparable or superior to P a for large n(n ^30) with P B be- 
ing somewhat superior when n is small, say n < 10. 



The numerical integrations for this study were performed using the NBS software package 
DATAPLOT developed by Dr. J.J. Filliben, and the graphs were prepared using the same package. 
The authors wish to acknowledge the helpful suggestions for changes in the manuscript made by Dr, 
P. Smith and Mrs. M. Natrella. 
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8. Appendix A 

The approximation P a given in (3.1.9) involves the computation of L[a,b,Q) defined as 

CO CO 

Ha,b,g) = J" jg(z,y,Q)dydz. 

The computation of Ha,b,g) can be reduced to a single variable integration. When a and b are both 
positive [18], 

Ha,b,g) = w- J exp {— Ua 2 +b 2 —2ab cos w) cosec 2 w)<Jw 

£" arc cos q •» 

The following recursion relations hold: 

Li-a,b, Q ) = -Ua,b-Q) + \[l-h{b)] 
L(arb,Q) = -L{a,b,- Q ) + \ [l-A(a)] 
Li-a -b, e ) = Ua,b, e ) + | [hia)+k[b)] 
where AW = J exp{ - t 2 /2 )dt. 
The approximation P a can be computed for all values of a,b and Q using the foregoing equations. 
P a = W-a)-L[a,b,Q),a>0,b>0 
P a = <M-a) -*(-*>) + H-a,b,-Q),a < 0, b > 
P H = Ua,-b-Q}, a >0,b<0 
P a = Wb)-Li~a,-b, Q ),a<0,b<0 
where *(*) = J" exp(-f 2 /2|d*. 

.OO 

9. Appendix B 

Asymptotic independence of the sample mean and the (n _ /c) ln extreme statistic. 

Let Xj,*", X n be i.Ld with a p.d.f. /|x|. Denote the c.d.f. of the X's by Fix). Assume that X's have a 
finite mean pi and finite variance 2 . Let Xmi < ••• < X( n j be the order statistics. 
The conditional density of X^, m ", X (n> given that X (n _ fc | = x^_^ is given by 
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n— k — 1 b 

(n-fe-ll! n /(x (/) ) fc! n /(*(,-,) 

L, = 1 n-ft+1 ... 

'""*' n-ft-1 fc "» 

[ffe ( „. A ,)J U-fUteHt,)! 

Clearly, given that X( n _^j = «j n _^i, the joint conditional density may be regarded as the joint density 
of two dependent samples [Fj,*", F n _^_j} and (Jfj,"*, fP^}, where the F-sample has a p.d.f. 



^In-A)' 
= , iIx>x H , 



(2) 



and the ^-sample has a p.d.f. 



(3) 



(4) 



fix) 
# M = JZFT1 7 ' ltx>x in-k) 

= , Ux<x {n _ k) 

THEOREM. For every fixed k, \^i(X— p) is asymptotically independent of X| n _| c i as n -*■ ™. 
PROOF: Rewrite X m terms of the K's and the Ws. We obtain 

yfiJX -f) _ yfe-A-1 JY-ji) s^j^l (W-ybk Xjn-krH 

From (2) we have 

JxdFix) JdFbc) (5) 

Making use of (4) and (5), and letting A be the value of EY f . with X^.^j replaced by X( n _£), we get 

yfalX -/v ) vSi-*-l [F-EYj) yfe-fc-1 



| (n-A-1) «A-fi) , (y ~^ ^_ + X \n-W 

Since k is fixed, clearly [W— \i) A/a\/ n_ *0 in probability as n -*■ «>. To prove the theorem we need the 
following two lemmas in which we show that the second and the fourth terms tend to zero in probabili- 
ty. Then the theorem follows from the fact that the first term converges in distribution to iV(0,l > which 
is the "unconditional" limiting distribution of \fii[X ~y). 
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Lemma Usn-«, 



X,I - k) -OfaP. 



PROOF: For every £ > and for a fixed A, it follows from the Chebyehev inequality that 



vfi 



>£ 1< E(X {n _ ki )2 l £<ma*X.2) 

|^ 5 ^ ne. K]<n 



Let Yj = X ; 2 and H{y) = P\Yj > y]. 
Following a proof in Chung (1960), 



Pfmax Yj > y] = 1 - [tf(y)] n > n[l-H\y)] 
Kj<n 



and 



^£fmaxX ; .2| = I j {i-[Hiy)]"} dy > J [l-ff(y)]<fy < °° 
\<j<n " 

On the other hand, 

- E\maxXf} = f J u n-i du dy. 
U Kj<n H{y) 

Sinee the expectation is finite, we can take the limit asn - *" under the integral sign. As a result 

1 °° 1 

]im ± E [max X ; 2 ] = / / Urn u"" 1 du dy = 

n-*<= Kj<n H{y) n-*« 

LEMMA 2. For a fixed k, « k < n-1 , 



00 

*i-k-l / x dF{x) -OinPasn-*". 

X (n-k» 

PROOF: Since 

00 

X (n-fct 

oo 

+ ^k^ J [1-FM] dx, (8) 

X (n-fc) 
we will show that each term on the right side of (81 converges in probability to zero. 
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Set 



z= 1 — - (9) 

n—k—l 



Then 



P{i n -k-l)[l-F(X {n _ k) )]>x} = J. (") V(l-z)»-' -e"* if 



(10) 



i=n— A 



We see thatX( n _^)/\/n— k— I -*■ in Pas shown in Lemma 1 and («— A— 1) [l-f(X( n _^))] converges in 
distribution as shown in ( 1 0). Thus, the first term on the right side of (8) tends to zero in P. 
Finally, to show that the last term in (8) tends to zero in P, write this term as 



vfe-A-1 J [l-F{X)]dx = vf B -ft-lKl-F(X faHfc) ) U([l-F\x)Wx)Wl-FVC 

X (n-k) I X { n -k) 



ln-ft)' | 



Clearly, the part in braekets tends to zero in P ean be seen by the application of the L'Hospitals's 
rule to it. 
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